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ABSTRACT 



Context. The X-ray spectra of Active Galactic Nuclei (AGN) carry the signatures of the emission from the central region, close to the 
Super Massive Black Hole (SMBH). For this reason, the study of deep X-ray spectra is a powerful instrument to investigate the origin 
of their emission. 

Aims. The emission line most often observed in the X-ray spectra of AGN is Fe K ff . It is known that it can be broadened and deformed 
by relativistic effects if emitted close enough to the central SMBH. In recent statistical studies of the X -ray spectra of AGN samples, 
it is found that a narrow Fe line is ubiquitous, while whether the broad features are as common is still uncertain. We present here the 
results of an investigation on the characteristics of the Fe line in the average X-ray spectra of AGN in deep Chandra fields. 
Methods. The average spectrum of the AGN is computed using Chandra spectra with more than 200 net counts from the AEGIS, 
Chandra Deep Field North (CDF-N) and Chandra Deep Field South (CDF-S) surveys. The sample spans a broader range of X-ray 
luminosities than other samples studied with stacking methods up to z ~ 3.5. We analyze the average spectra of this sample using 
our own averaging method, checking the results against extensive simulations. Subsamples defined in terms of column density of the 
local absorber, redshift and luminosity are also investigated. 

Results. We found a very significant Fe line with a narrow profile in all our samples and in almost all the subsamples that we 
constructed. The equivalent width of the narrow line estimated in the average spectrum of the full sample is 74 eV. The broad 
component of the Fe line is significantly detected in the subsample of AGN with L < 1.43 x 10 44 erg s and z < 0.76, with an 
equivalent width 108 eV. 

Conclusions. We concluded that the narrow Fe line is an ubiquitous feature of the X-ray spectra of AGN up to (z) ~ 3.5. The broad 
component of the Fe line is significantly detected in the X-ray spectra of the AGN with low luminosity and low reshift. 

Key words. Galaxies: Active -X-rays: galaxies 



1. Introduction 

The X-ray spectra of AGN are emitted from the innermost 
regions of the central engine, close to the central SMBH. 
Accordingly to the currently accepted model described in detail 
by |Shakura et al. (19 74) AGN are powered by accretion to the 
SMBH, with large amounts of potential energy released in the 
accretion disc. Assuming an optically thick, geometrically thin 
. accretion disc, the optical-UV primary emission of the AGN is 
$_i ' explained as the result of the thermal emission from the accre- 
tion disc. 

It is well known that the emission in hard X-rays is domi- 
nated by a powerlaw with F ~ 2, and the most widely accepted 
interpretation for it invokes a hot plasma surrounding the accre- 
tion disc that makes inverse Compton scattering of the thermal 
optical-UV photons from the accretion disc (Ha ardtetal. 199 l"b . 

Part of the primary emission is reflected by the accretion 
disc producing Compton reflection (especially significant above 
10 keV) and several fluorescence lines, the most important one 
being the Fe K a line at 6.4 keV for neutral Fe (and 6.7-6.9 
keV for Fe in the most highly ionized states), as described in 
IMattetal. (1991)1 and |George & Fabian (1991 )| Reflection can 
also occur at larger scales, in the torus. The line profile gives 



valuable information about the emitting processes and the re- 
gions where they occur. 

The study of the trend of the narrow Fe line equivalent width 
(EW) with continuum luminosity plays a fundamental role in 
the investigation of the link between the continuum emitting 
region and the reflecting region. A decrease of Fe line inten- 
sity with increasing continuum luminosity was found in Ginga 
observations of AGN ( Iwas awaetal. 1993b . The trend (the so 
called "Iwasawa Taniguchi effect", hereafter "IT effect") was 
widely confirmed for both broad and narrow Fe lines by ob- 
servations of ASCA, by Nandra et al. (1997 )] and more rec ently 
for the narrow Fe lines, for example in Bia nchi et al. (2007)| and 
|Page et al. (2004)| using XMM-Newton data. Although the IT ef- 
fect is widely proved, the reason for the anti-correlation between 
the line and the continuum intensity is still under investigation: 
|Nandra et al. (1997)| attributed the IT effect to the ionization that 
can be important in high luminosity. 

As the torus and the disc can both contribute in different mea- 
sure to the observed Fe line, the resulting profile will be com- 
plex, and in particular it will depend on the inclination angle of 
the disc with respect to the line of sight. The estimated intrinsic 
column density in X-rays (Nh) is, in this sense, connected to the 
line properties: a increase of the Fe line intensity with increas- 
ing Nh was found in a sample of AGN observed by Ginga and 
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Asca (Gil li et al. 1999b . In that work, the Fe line in the spectra of 
the absorbed AGN (type 2 AGN) was more intense than the line 
found in the unabsorbed ones (type 1 AGN). This was predicted 
by Ghisell ini et al. (1994)| who performed a theoretical analysis 
of the Fe line in a configuration with torus and accretion disc 
on the same plane. In this picture, the line intensity produced by 
reflection and transmission in the torus increases with the Nh, 
while the continuum emission is depleted, leading to an increase 
of the line equivalent width for the highest Nh- 

It should be noted that the lines produced by reflection in the 
torus should have a narrow and approximately Gaussian profile, 
while the lines produced in the accretion disc should be broad- 
ened and deformed by the effect of the newtonian movement 
of the gas in the accretion disc around the SMBH. It is also 
possible to detect in the lines the relativistic features, such as 
a red wing due to gravitational redshift, among other features, 
explained in detail in Fabian ^t al. (2000)1 These features come 
from the innermost regions of the accretion disc, as the relativis- 
tic effects are due to the strong gravitational field close to the 
central SMBH. From the study of these features it is possible to 
obtain the inner radius of the accretion disc, which is expected 
to be 6 gravitational radii (R g ) for a non-rotating BH, called a 
Schwarzschild BH, d own to 1.23 R g for a maximally rotating 
BH, called a Kerr BH dBardeen at al. 1972K 

Early results from ASCA showed broad and deformed lines: 
|Nandra et al. (1997)| performed a spectral analysis of a sample of 
local Seyferts, finding broad components in about 65% of their 
sample. Unambiguous broad lines were found also in more re- 
cent well exposed XMM-Newton observations of AGN (e.g. in 
INandra et al. 20061 iBraito et al. 20071 and IFabian et al. 20021 

It has been proved by Guai nazzi et al. (2006) that X-ray 
spectra with good statistical quality are necessary to detect any 
broad Fe line component, which would otherwise be hidden un- 
der the noise. In this context statistical methods, like summing 
(stacking) or averaging the spectra, have been recently intro- 
duced in X-ray astronomy to allow the study of large AGN sam- 
ples including the low quality spectra that otherwise could not 
be well analyzed individually. 



The first of such works is that of Streblyanska et al. (2005) 



who performed a stacking analysis on the XMM-Newton obser- 
vation of the Lockmann Hole, finding a large line EW and broad 
line profile. 

|Brusa et al. (2005)| stacked the spectra of CDF-N (2 Ms ex- 
posure time) and CDF-S (1 Ms exposure time) surveys by the 
Chandra satellite. They computed the stacked spectra in bins 
of redshift aiming at characterizing the Fe line emission of the 
sources in the X-ray background up to z ~ 4. They found an in- 
tense and apparently broad 6.4 keV Fe line with an EW consis- 
tent with the results of Streblyanska et al. (2005) To explain the 
red component, they argued that it can strongly depend on the 
modeling of the underlying continuum and that a spurious red 
wing might be produced by the contribution of absorbed spectra 
at different redshifts. 

In recent studies of X-ray spectra using this approach, find- 
ing broad lines with relativistic profiles has been proved to be not 
very common, for example in Guai nazzi et al. (20 06) in the anal- 
ysis of AGN spectra observed with XMM-Newton. They found 
that 25% of their sample had relativistic lines; this percentage is 
50% for higher signal-to-noise ratio, selecting only spectra with 
large numbers of counts. Moreover, they found the strongest rel- 
ativistic profiles in low luminosity objects. 

|Corral et al. (2008)| average d the XMM-Newto n spectra of 
type 1 AGN, using the AXIS dMateos et al. 20051 1 and XWAS 
(Mate os et al. 20 10) samples, up to redshift ~ 3.5. They found 



a narrow Fe line, while no clear evidence of a broad line was 
found. The reason for the discrepancies between this work and 
that of Streblyans ka et al. (2005)| can be explained consider- 
ing the differences of the samples, as the |Corral et al. (2 008) 
sample has higher luminosities (and therefore lower predicted 
EWs, following the IT effect) and less spectral counts, in- 
cluding more noisy sources. Another reason for the discrep- 
ancy can be found in the method, since [Corral et al. (2 008) 
estimate the continuum shape using simulations, while in 



Streblyanska et al. (2005) it was not constructed and subtracted 



in the same way and this could introduce some uncertainties. 



Moreover, Streblyanska et al. (2005)1 binned the spectra before 



stacking them, and |Yaqoob (2006)| showed that this procedure 
can introduce features like a broad red tail in an emission line. 

A stacking analysis of a deep and complete sample of 507 
AGN with z < 4.5 defined from the 2XMM catalogue was per- 



formed by Chaudha ry et al. (2010)| They were able to character- 
ize the properties of the stacked spectrum of the AGN, such as 
the Fe line shape and the dependence of its intensity with X-ray 
continuum luminosity and the redshift. They found a clear ev- 
idence of a narrow neutral Fe line, and they confirmed the IT 
effect in AGN over a broad range of redshift. 

Recently, a stacking analysis on XMM-Newton X- 
ray spectra of the COSMOS sample was performed 
dlwasawa et al. 20PT1) : they found an excess on the high- 
energy side of the Fe line, interpreted as the convolution of 
narrow lines from ionized Fe. 

In summary, according to recent results, while the nar- 
row Fe lines are commonly detected in the spectra, clear ev- 
idence for relativistic lines is rare. The paucity of relativis- 
tic fines is a problem, as they are expected in the accre- 
tion disc scenario. There are a number of possible solutions 
for this problem. The first one invokes the presence of ion- 
ized discs (Iwasa wa et al. 20PT1 IMatt et al. 19961 ); the second 
solution explains it with the presence of fastly spinning BH 
dlwasawa et al. 19961 and IFabian e t al. 2002): in this case the 
lines would be very broad and sometimes difficult to separate 
from the continuum. The last explanation for the lack of broad 
lines postulates truncated discs, as found by Matt et al. (2005) 

We present in this paper a stacking analysis of a com- 
prehensive sample of absorbed and unabsorbed AGN in the 
deepest Chandra surveys: AEGIS and Chandra Deep Fields. 
We followed the method presented by |Corral et al. (2008)] and 
we accurately tested and adapted it for Chandra spectra, in 
order to carefully check our results. The Chandra sample in 
this work allows us to explore the properties of AGN with 
the best statistics ever reached at high redshift, as it is deeper 
than the XMM-Newton samples that have previously been an- 
alyzed with stacking techniques, such as |Corral et al. ( 2008) 
and Chaudhar y et al. (2010)| (see Sect. 2.2). In comparison with 
the previous stacking analysis of the Chandra samples in 
Brusa et al. (2005), we include the sources from the AEGIS-X 
survey, improving the statistics of high luminosity AGN, and we 
use the more recent (2 Ms) observation of the CDFS (Sect. 2.1). 

The paper is organized as follows: the properties of the sam- 
ple are described in Sect. 2; the method in Sect. 3; discussion of 
the results is in Sect. 4; we summarize the conclusions in Sect. 


Across this paper, we adopt the cosmological parameters: 
H o=10 km s _1 Mpc' 1 : Q^0.7 dKomatsuet al. 201 II) . 

In this paper, all the counts refer to the net number of counts 
between 2 and 12 keV rest-frame; all the luminosities are cal- 
culated between 2 and 10 keV rest-frame and are corrected 
for Galactic and intrinsic absorption using the fits to the in- 
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dividual spectra (see below). We have used XSPEC v. 12.5 
( [Arnaud, 19 96) for spectral analysis. 



2. X-ray samples 

The X-ray results presented in this paper are obtained from the 
Chandra observations of the AEGIS, CDF-N and CDF-S sur- 
veys. 

2. 1 . Parent Chandra surveys 

The AEGIS-X survey is a Chandra survey of the Extended 
Groth Strip (EGS) region, designed primarily for studying the 
co-evolution of black holes and their host galaxies. The data 
we used are the result of 8 contiguous Advanced CCD Imaging 
Spectrometer (ACIS-I) pointings, each with exposure 200 ks, to- 
talling 1.6 Ms. The survey covers a total area of approximately 
2302 arcmin 2 a in a strip of length 2 degrees. The total number of 
identified sources with spectroscopic redshifts is 409. The data 
reduction and point source detection algorithms used to analyze 
these data are described in |Laird eta l. (2009) 

The Chandra Deep Field-North (CDF-N) survey is one of 
the deepest (2Ms with ACIS-I) 0.5-8.0 keV surveys ever made: 
nearly 600 X-ray sources are detected over 448 arcmin 2 . The 
total number of sources identified with spectroscopic redshifts is 
307. Details of data reduction and the point source catalog are 
described in the paper by Alexande r et al. (2003)[ 

The Chandra Deep Field South (CDF-S) is the deepest and 
most sensitive observation obtained with Chandra. We used 
the 2 Ms survey, covering an area of 436 arcmin 2 . Several 
hundred X-ray point sources were detected. The number of 
sources identified with spectroscopic redshifts is 152. A de- 
tailed description of the survey can be found in the paper 
|Luo et al7 (2008) It has been recently extended to 4 Ms, as de- 
scribed in|Alexander et al. (201 1)[ and|Xue etal. (201 1)| 
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Fig. 1. Distribution of the net counts between 2 and 12 keV 
rest-frame. Filled histogram: CDF-S; histogram with horizontal 
stripes: AEGIS, histogram with vertical stripes: CDF-N. 
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Fig. 2. Redshift distribution. Filled histogram: CDF-S; his- 
togram with horizontal stripes: AEGIS, histogram with vertical 
stripes: CDF-N. 



2.2. Sample definition and properties 

Our main purpose is to study the average spectrum of AGN in the 
2-12 keV energy band (rest- frame), and to determine the shape 
of the Fe line at ~6.4 keV. The X-ray spectra of our sample have 
few counts and do not allow a detailed analysis of the individ- 
ual spectra, except for a handful of sources. For the detection 
of any relativistically broadened component, as underlined by 
|Guainazzi et al. (2006)| it is necessary to have good quality spec- 
tra, hence we need to improve the signal-to-noise ratio. 

To increase the probability of detecting any broad compo- 
nent that otherwise would be hidden by the noise, we included 
in our analysis only those sources with more than 200 counts. 
We discuss in Appendix A the results corresponding to lower 
thresholds (50 and 100 counts). We further selected only sources 
with spectroscopic redshifts. Finally, we excluded from our anal- 
ysis also two sources with more than 10000 counts individu- 
ally: CDFS-056 (RA: 53.112, DEC: -27.685), CDFN_141 (RA: 
189.1, DEC: 62.383). These sources have such a strong signal 
in their spectra that they would largely dominate the average 
stracked spectrum. 

Our full sample contains 123 sources with spectroscopic red- 
shifts and more than 200 counts each (but less than 10000 cts), 
with 70667 counts in total (see Table [TJ- The distribution of the 
net counts per source for this full sample is shown in FigQ] The 
redshift distribution of the sources is shown in Fig. [31 where we 



can see that our sample spans a broad range of redshifts from 
z = up to z = 3.5. 

Source and background spectra, and their ancillary 
and response matrices were extracted for each source 



(Dig by-North, 201 1[>, using the tool AC IS -EXTRACT, version 
2008-03-04, |Broos et al. (2010)1 Before spectral fitting, all the 
spectra were grouped with a minimum of 5 counts per bin, 
and we applied the modified Cash statistics discussed in the 
XSPEC pagesQ. This grouping was adopted because we used the 
same method for all the sources with more than 50 counts (see 
Appendix [A), for which the more common x 2 statistics would 
not be adequate. 

It should be noted that in this step we did not aim at making 
a detailed analysis of the individual spectra but rather only at ob- 
taining a set of parameters that we use to unfold the spectra for 
all the instrumental effects (see Sect. [3}. The model is a single 
power law modified by Galactic absorption (with fixed values at 
1.4 xl0 20 cirr 2 for AEGIS, 0.772 xl0 20 cirr 2 for CDF-S, 0.993 
xl0 20 crrT 2 for CDF-N) plus intrinsic absorption at the redshift 
of each source. We left the slope of the power law, its normal- 
ization and the intrinsic column density as free parameters. The 
fit was performed in the rest-frame energy range between 1 and 
12 keV, minimizing the contribution from any putative soft ex- 
cess present in the sources. The resulting distribution of intrinsic 



http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/wstat.ps 
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Fig. 3. Distribution of intrinsic column density in units of 10 
cm -2 . Filled histogram: CDF-S; histogram with horizontal 
stripes: AEGIS, histogram with vertical stripes: CDF-N. We 
grouped the spectra with log(NH) < 10 20 with the spectra having 
log(N H ) = 10 20 . 
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column densities is shown in Fig|3] Most of the highly absorbed 
sources belong to CDF-S, as indicated by the higher average col- 
umn density for that sample in Table Q] 

2.3. Subsample definition 

We show in Fig. [4] the distribution of the X-ray luminosities of 
the sources corrected for Galactic and intrinsic absorption, cal- 
culated from the fit to the indvidual spectra. The average lumi- 
nosity of the AEGIS sources is one order of magnitude higher 
than the average luminosity of CDF-N and CDF-S sources 
(TableQ). 

We compared the distribution of our ful sample in the 



luminosity-redshift plane with that of Streblyanska et al. (2005) 
in the Lockmann Hole (Fig. [5): we sample better that plane at 
Z < 3 for all luminosities. Compared to the |Corral et al. (2 008) 
sample (XMS-XWAS) (Fig.|6]l we cover about one order of mag- 
nitude lower luminosities for 0.5 < z < 3. Finally, with respect 
to Brusa et al. (2005) the addition of the AEGIS sample and the 
restriction to sources with high numbers of counts, allows us to 
span a broader range of luminosities at similar z. In summary, 
our sample combines a unique combination of deep coverage 




redshift 



Fig. 5. Comparison of our full sample (black stars) with the deep 



XMM-Newton sample of Streblyanska et al. (2005) (LH, red tri 



angles), in the luminosity-redshift plane. 




Fig. 6. Comparison of our full sample (black stars) with 
the medium XMM-Newton sample of |Corral et al. (2 008) 
(XMS+XWAS, red triangles), in the luminosity-redfhit 
plane. For clarity, we omitted one XMS+XWAS source of 

erg s _1 . 



XMS+XWAS at z = 2.34, L - 1.28 x 10 47 ■ 



with high z sources, a broad span of luminosities and good statis- 
tics. 

One of the main aims of our work is to understand the de- 
pendence of the spectral properties of the AGN on luminosity, 
redshift, and intrinsic absorption. For this reason, we performed 
the analysis of the stacked spectra not only for the full sample, 
but also for subsamples defined in terms of column density, lu- 
minosity, and redshift. Except for Nh (see below) the bins have 
been designed to have a similar numbers of total counts each. 
The characteristics of each subsample can be seen in TableQ] 

- Intrinsic A^: The threshold column density was set at 
log(A r n/cm _2 ) = 21.5 as the absorption features of the X-ray 
spectra are detected in the AGN with log(A r n/cm _2 ) > 21.5. 
Local absorbers with log(Nn/cnT 2 ) = 21.5 are opaque for 
E < 1 keV and transparent for E > 1 keV. The transition 
energy between the transparent and opaque regime linearly 
grows with jVh, until A^h ~ 10 24 cirT 2 , where the Compton 
Thick regime dominates. We show in Fig.|7]the distribution 
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Fig. 7. Distribution of the unabsorbed (red triangles) and ab- 
sorbed (black diamonds) subsamples (see Sect. 2.3) in the 
luminosity-redshift plane. 
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Fig. 8. Distribution of the L — z subsamples (see text) in the 
luminosity-redshift plane. The dividing point is z — 0.76, L = 
1.435 x 10 44 erg s _1 (marked by the horizontal and vertical solid 
lines). The three blue stars mark the sources that allow defining 
distinct L - z subsamples (see text). 
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Fig. 9. Average observed spectrum of the full sample fitted with 
a powerlaw between 2 and 5 keV. 



results are very similar and we have used the latter definition 
of distinct areas as our default. 



3. Analysis method 

3. 1 . Averaging method 

Our averaging method was presented for the first time in 
ICorral et al. (2008)| for the XMM-Newton spectra of XMS- 
XWAS surveys; it was then also used for the study of the XBS 
sample by |Corral et al. (201 1)1 We adapted and tested it for the 
Chandra spectra. 

After fitting the spectra as described in Sect.2.2, we read into 
XSPEC again each un-binned, background-subtracted spectrum, 
we applied the corresponding best fit model and we extracted 
the unfolded spectra taking into account the detector response 
(euf spec command in XSPEC). In this step it is possible that 
some distortions to narrow spectral features were introduced, 
such as to the Fe emission line. In order to quantify this effect, 
we performed extensive simulations of the continuum and of the 
Fe Line as described in Sect. 3. 2 and 3.3. Once we obtained all 
the unfolded spectra, the following step was to apply the correc- 
tions described in |Corral et al. (200 8) that we can summarize as 
follows: 



of the absorbed and unabsorbed subsamples just defined: the 
two subsamples have a similar distribution in the luminosity- 
redshift plane, so they are fair representations of the AGN 
with and without absorption, with no other parameters play- 
ing an important role (see also Table[TJ. 

- Luminosity: We separated high and low luminosity sources 
using a threshold L = 8 x 10 43 erg s _1 

- z: We built a low-z and a high-z subsample with a threshold 
z= 1.005 

- L—z: Initially we divided the full sample in three subsamples: 
(i) z < 0.76 (hereafter low z-low L) (ii) z > 0.76 and L < 
14.35 x 10 43 erg/s (hereafter high z-low L) (iii) z > 0.76 and 
L > 14.35 x 10 43 erg/s (hereafter high z-high L). In the low z- 
low L subsample there are only 3 sources above L < 14.35 x 
10 43 erg/s (with ~4000 total counts), so we have also used a 
second version of this subsample using only sources with z < 
0.76, L < 14.35 x 10 43 erg/s so that this subsample and the 
two high z subsamples cover distinct z and L intervals. All 



1 . Correction for Galactic absorption 

2. Shift the spectra to a common redshift frame 

3. Re-normalize the spectra with respect to the integrated flux 
between 2 and 5 keV restframe 

Once we obtained the de-absorbed and re-normalized spec- 
tra at rest-frame, we re-binned each spectrum. For simplicity, we 
used the same binning for all the average spectra, applying the 
bin definitions of the absorbed sample: this binning was con- 
structed in order to have at least 1000 net counts in each bin of 
the co-added spectrum. In order to maximize our ability to de- 
tect a narrow Fe line, we always centered one bin at 6.4 keV. 
We finally averaged the spectra using the unweighted arithmetic 
average. 

We show in Fig.|9]the average spectrum of the total sample. 
If we fit the 2-5 keV range with a powerlaw and extrapolate it 
to the full range, an excess can be clearly seen around the ex- 
pected position of the Fe emission line. However, any features 
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Fig. 10. Average observed spectrum of the full sample (data 
points) with the average simulated continuum (continuous line) 
and the one sigma confidence limits (dashed lines). 



found in the average spectra must be carefully evaluated, as the 
averaging procedure itself may introduce some distortions. For 
this reason, we computed the continuum using simulations, as 
described below. 

3.2. Simulations of the continuum 



Previous works, as for example Yaqoob (2007) and 
|Corral et al. (2008)] have shown that the unfolding and the 
averaging process can distort the shape of the spectrum. 
Therefore, it is important to take this into account before 
deriving conclusions from the spectral analysis of the average 
spectrum. 

In order to characterize the underlying continuum of our 
sample, we made 100 simulations of each source using the best 
fit parameters of the continuum model (see Sect. 2.2). To each 
of these 100 simulated samples we applied the same method as 
the one used for the observed sample (spectral fitting, correcting 
for response, correcting for Galactic absorption and z, rescaling 
and averaging). After this, we represented our continuum with 
the median of the 100 averaged simulated continua. We decided 
to use the median and not the arithmetic average as it is a more 
robust estimate of a central value and in particular it is less sen- 
sitive to extreme values. Hereafter we will call this the simulated 
continuum. 

We represent in Figs.[10][TTJ[T2][T3] and[T4lthe average ob- 
served spectra (dots with error bars), along with the simulated 
continua (solid lines), and the 1 sigma confidence limits (dashed 
lines). The latter correspond to the 16th and 84th element of 
the sorted fluxes of the simulated spectra at each bin. Features 
around the expected location of the Fe K line are conspicuous in 
all the subsamples. 

Before attempting to draw any conclusions from those fea- 
tures, we have characterized the effect of our averaging method 
on narrow lines using simulations, as described in the next sec- 
tion. 

3.3. Simulations of the Fe Line 

We performed simulations of unresolved Fe lines in order to 
study how the spectral resolution of ACIS-I and the averaging 
process widen the narrow spectral features in X-ray spectra. 

We simulated high signal to noise spectra (one simulation 
for each spectrum, without Poissonian noise) using a powerlaw 



with T fixed at 1.9, with unit normalization, and a zero sigma 
Gaussian. The input EW was 200 eV and energy fixed at 6.4 keV, 
at the z of each source. We made the correction for the response 
(using a power law with Gamma fixed at 1.9) and we finally 
computed the average spectrum with the same treatment applied 
to the real spectra. 

We show in columns r s , m /,„ c , cr^ m u ne and EW of Table Q] the 
results of fitting a powerlaw plus a Gaussian to these simulations 
for some of the samples. The EW and powerlaw slope are recov- 
ered very well. The width of the line has become cr ~ 120 eV. 
Therefore, our resulting "instrumental resolution" is around 120 
eV, i.e., any real detected feature around rest-frame 6.4 keV 
should be wider than this value and its actual width will be the 
convolution of its intrinsic width with this value, roughly added 
in quadrature (see Eq. 1 below). 

We show in the bottom-left corner of Fig. [14] the expected 
profile of an unresolved 6.4 keV line added to the simulated 
continuum for the corresponding subsample. In general, the fea- 
tures we observe around that energy appear to be compatible 
with such an unresolved line, except perhaps for the low z-low 
L, high z-high L and high L subsamples. We will quantify these 
qualitative impressions in Section l4~2l 



4. Results 

In this section we study the statistical significance of the Fe line 
and its properties in several steps. We first developed a model- 
independent way to calculate the significance of the line detec- 
tion (Sect. 4.1), then we studied again its significance and char- 
acteristics fitting the average spectra using a single Gaussian to 
fit the Fe line (Sect. 4.2). In Sect. 4.3 we study the dependence 
of the line on z and X-ray luminosity. Finally, in Sect. 4.4 we try 
complex line models with those subsamples that show signifi- 
cant broad line profiles. 

4.1. Significance from the simulations 

In order to estimate the significance of the excess observed in the 
Fe K line region, we calculated the percentage of average sim- 
ulated spectra (from the simulations of the continuum described 
in Sect. 3.2) having a flux in the Fe K region which is lower 
than the flux of the observed spectrum. We should note that this 
method does not depend on a modeling of the continuum (as in 
the calculation of the equivalent width, in the next Section). The 
Fe line flux was calculated in three regions centered at 6.4 keV, 
with half- widths of 0. 1 keV, 0.2 keV, and 0.4 ke V respectively. 
The significances of the Fe line estimated in this way for the full 
sample and all subsamples are in Table [2] 

The narrowest interval was chosen at the limit of the instru- 
mental broadening, a significant excess found just in this range 
and not in the other two would in principle correspond to a detec- 
tion of a narrow line. On the other hand, increasing the width in- 
cludes more flux from a putative broad line, but it also increases 
the noise, hence a broad component would be more difficult to 
detect. We use the +0.2 keV wide option (column S 2 in Table 
|2]l as our fiducial value, since it maximizes the signal-to-noise 
ratio for unresolved and moderately broad lines. We can see that 
the line is strongly significant (> 98 %) in the full sample and in 
almost all the subsamples, with the exception of the high z-high 
L one. In general, the significance is higher for low z and low L. 

In all cases, the significance stays about the same or de- 
creases when AE increases. This is either an effect of the ex- 
pected lower signal-to-noise ratio for wider intervals, or might 
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log (N H cm- 2 ) > 21 .5 log (N H cm- 2 ) < 21 .5 




2 5 10 2 5 10 

Energy (keV) Energy (keV) 



Fig. 11. Average observed spectra (data points) of the absorbed (left) and unabsorbed (right) subsample with their average simulated 
continua (continuous line) and the one sigma confidence limits (dashed lines). 



z > 1 .005 z < 1 .005 




2 5 10 2 5 10 

Energy (keV) Energy (keV) 



Fig. 12. Average observed spectra (data points) of the high-z (left) and low-z (right) subsamples with their average simulated continua 
(continuous line) and the one sigma confidence limits (dashed lines). 



L < 8 1 43 erg/s L > 8 1 13 erg/s 




Energy (keV) Energy (keV) 



Fig. 13. Average observed spectra (data points) of the high-L (right) and low-L (left) subsamples with their average simulated 
continua (continuous line) and the one sigma confidence limits (dashed lines). 



indicate the presence of a wider profile in the cases where the 
significance does not decrease noticeably. 



4.2. Spectral fits on the full sample and subsamples 

After having made the above mentioned calculation of the sig- 
nificance of the Fe line, we performed the spectral analysis. To 
do this, we represented the continuum using the simulated con- 



tinuum as a table model in XSPEC. During the fits described 
below, we have aways left its normalization free to vary, and its 
value returned in XSPEC is always around one. We added to the 
table model the required components when we found residuals. 
We used the goodness of fit criterion and the confidence inter- 
val for a given parameter corresponding to a Ax 2 = 2.7 1 (90 % 
probability). The fit results are shown in Table|5] 
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L > 1 .43 1 " erg/s, z > 0.76 




Energy (keV) Energy (keV) 



Fig. 14. Average observed spectra (data points) of the L — z subsamples with their average simulated continua (continuous line) and 
the one sigma confidence limits (dashed lines). For comparison, in the bottom-left panel we represent with a dotted line the result 
of adding to the simulated continuum a Gaussian centered at 6.4 keV, with <x=120 eV, as obtained from the simulations of the Fe 
line (Sect. 3.3). 



Table 2. Significance of the Fe line estimated from the simula- 
tions for the full sample and the subsamples. 



Sample 


Si 


s 2 


S3 




(1) 


(2) 


(3) 


Total 


100 


100 


99 


log(N H ) > 21.5 


99 


99 


98 


log(N H ) < 21.5 


100 


100 


100 


L 43 < 8 


100 


100 


100 


L 43 > 8 


97 


98 


96 


z < 1.005 


100 


100 


100 


z > 1.005 


98 


98 


98 


L 43 < 14.3, z < 0.76 


100 


100 


100 


z < 0.76 


100 


100 


100 


L 43 < 14.3, z > 0.76 


99 


98 


91 


L 43 > 14.3, z > 0.76 


88 


88 


86 



Notes. L 43 : luminosity in units of 10 43 erg s -1 . Si, S2, S3: percentage 
of average simulated spectra (from the simulations of the continuum, 
see Sect. 3.2) with the line flux lower than the observed spectrum in the 
intervals: 6.3-6.5 keV (1); 6.2-6.6 keV (2); 6.0-6.8 keV (3). 

The cr we mention in the remainder of this work (including 
Table [5]), refers to the intrinsic width of the line, taking into ac- 
count the broadening introduced by our method 



where cr b s is the total width returned by XSPEC and cr sim is the 
instrumental width obtained from the simulations of an unre- 



solved Fe K line (120 eV, see Sect. 3.3.). During the fits, we 
forced the width of the Gaussian line to satisfy the condition: 
cr bs > o~sim (where we used the average value of cr sim reported 
in TableQ] that is approximately 120 eV). 
We made the fits in several steps: 

- Fixed o"=0 and fixed centroid energy at 6.4 keV: this allows 
us to estimate the significance of the narrow component of a 
neutral Fe K line (first line for each sample in Table 3) 

- Fixed centroid energy at 6.4 keV and free <x: studies the sig- 
nificance of a possible broad component and constrains its 
width (second line for each sample) 

- Free centroid energy and fixed <x=0: considers the presence 
of an ionized narrow Fe component and estimates its centroid 
energy (third line for each sample) 

- Free centroid energy and free cr: leaves all options open 
(fourth line for each sample) 

We calculated the significance of the Gaussian with the Ay 2 
corresponding to the fits with and without the line component, 
and we checked the corresponding probability using the incom- 
plete f3 function (according to lPress et al. 20071) . The result is in 
Col. 4 in Table 3. In the same way, we estimated the significance 
of allowing the width and/or the centroid energy to vary with 
respect to the baseline narrow and neutral line. 

We summarize here the results for our (sub)samples (see 
Tableland Figs. El OH El HI and[l9j : 

- Full sample: 

The spectrum of the full sample fitted with the table model 
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and the Gaussian with free energy and free a are shown in 
Fig US] We can see that the model fits the spectrum well in 
the Fe line region. The detected line is narrow. We do not 
observe any ionized Fe line component. 

- Absorbed and unabsorbed subsamples: 

The spectra and fitted models are shown in Fig [16] The con- 
tinuum of the unabsorbed sample shows F ~ 1.8, a common 
value for the AGN (see Table 1), and is flatter in the absorbed 
sample, as expected. The line is characterized by having a 
narrow profile for both subsamples. Only in the unabsorbed 
sample the centroid energy of the line might suggest a contri- 
bution from mildly ionized Fe. However, the centroid energy 
is consistent with 6.4 keV within the 90 % confidence level, 
and the significance of cr >0 from the Ay 2 is < 90 %. 

- High and low luminosity subsamples: 

The spectra and fitted models are shown in Fig. [17] The re- 
sults of the fits suggest a more intense Fe line in the low-L 
subsample. We discuss the dependence of the EW on red- 
shift and luminosity more in detail in the next Section. The 
detected line is narrow, and there is no significant ionized Fe 
contribution. 

- High and low redshift subsamples: 

The spectra and fitted models are shown in Fig. [18] The Fe 
line is more significant at low redshift than at high redshift 
and this can reflect the same trend just found with the lumi- 
nosity, since sources with higher L are preferentially found 
at higher z. The line profile is narrow in both cases. The fits 
with the centroid energy free shows the presence of a mildly 
ionized Fe contribution in the low redshift sample, but 6.4 
keV is within the 90 % confidence interval. 

- L-z subsamples: 

The strongest line significance is found at low z and low L, 
then decreasing with both redshift and luminosity. The line 
profile is broad only in one case, in the low redshift- low 
luminosity sample (>95.4 % probability). The line profile 
looks symmetric (see Figs. 14 and[T9l two left panels). This 
suggests reflection from regions in newtonian movement in 
the accretion disc, although we can not exclude the relativis- 
tic profile. Moreover, in the top-left panel there seems to be 
an excess in the narrow bin around 6.4 keV, which might 
come from an additional narrow component. We discuss in 
more detail the shape of this feature in Section I4~4l In the 
other two subsamples, the detected profile of the Fe line is 
narrow. We did not find any ionized Fe component. 

At high energies (~ 10 keV), we observed in some spectra an 
excess that can be interpreted as part of the Compton reflection. 
The paucity of counts at these energies does not allow us to make 
a more accurate assessment of this continuum component. An 
absorption feature was also detected, to some extent, in some 
of our average spectra around 7-8 keV: this can be explained 
as the Fe edge commonly found along the Compton reflection 
component. 

4.3. Dependence on redshift and luminosity 

As our results seem to suggest a dependence of the Fe line EW 
on redshift and on luminosity, for low z and low luminosity 
AGN, we checked this result in more detail to understand it bet- 
ter. The distribution of the line EW with the average z and the av- 
erage luminosity of the subsamples are shown in Fig. [20] In those 
figures, the same symbols refer to statistically independent sub- 
samples. On the contrary, subsamples with different symbols are 
not statistically independent and there will be significant over- 



Total sample 

0.1 | , , , 




Energy (keV) 



Fig. 15. Fit of the average spectrum of the full sample using the 
simulated continuum model + Gaussian (with free cr and free 
centroid energy) 



lapping among the sources included in each subsample. We have 
used the fits corresponding to an unresolved line fixed at 6.4 keV 
(first line for each subsample in Table 3). 

At first view the EW seems to decrease both with increasing 
z and X-ray luminosity. However, looking for example at the 
blue squares (corresponding to independent L-z subsamples), 
it is clear that a constant EW is consistent with the data points 
the upper limits and, hence, the trend, albeit suggestive, is not 
statistically significant. 

4.4. Analysis of broad lines 

We detected with strong significance narrow Fe K lines in all 
samples. Additionally, in the low L-low z subsample, the Fe line 
appears significanlty broadened with a significance above two 
sigma). We investigate here whether relativistic profiles are bet- 
ter fits than a simplistic broad Gaussian. 

We added to the simulated continuum the diskline model 
in XSPEC which describes a relativistically broadened emission 
line for an accretion disc around a Shwartzschild Black Hole. We 
fixed R ou t, the external radius of the accretion disc, to 1000 R g . 
We left R/„, the inner radius, free to vary, obtaining an upper 
limit for it: < 142 R g . As mentioned in Sect. 1, the inner ra- 
dius of the accretion disc is expected to be lower (6R g for a 
non-rotating black hole and 1.23 R g for a maximally -rotating 
one). The EW is consistent with the one obtained in the fit with 
the Gaussian (see Table 3): EW=152^[j eV. The significance of 
adding an Fe line with diskline to the continuum is > 99.99%, 
again consistent with the case of the Gaussian (Table 3). The fit 
gives ^ 2 /dof=20.53/29, to be compared with x 2 /dof= 19. 07/25 
for the Gaussian: the fit is not better with more free parame- 
ters. Summarizing, the diskline fits the Fe K line as well as the 
Gaussian: we do not significantly detect a definite relativistic 
profile in the line, although it cannot be excluded. 

Alternatively, the profile of the line could be characterized as 
the sum of both a broad and a narrow components since, in prin- 
ciple, we would expect a contribution from both a torus (narrow 
component) and disk (broad component) reflections. 

We have assessed the significance of this possibility by mak- 
ing a final test in the low L- low z sample: we added a broad 
Gaussian to the narrow Gaussian (i.e. combining models in the 
first and second line of Table 3). We fixed the centroid energies 
of the two Gaussians at 6.4 keV and we fixed the cr to zero in the 
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Fig. 16. Fit of the average spectrum of the absorbed (left) and unabsorbed (right) subsamples using the simulated continuum model 
+ Gaussian (with free cr and free centroid energy) 
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Fig. 17. Fit of the average spectrum of the low luminosity (left) and high luminosity (right) subsamples using the simulated contin- 
uum model + Gaussian (with free cr and free centroid energy) 
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Fig. 18. Fit of the average spectrum of the low redshift (left) and high redshift (right) subsamples using the simulated continuum 
model + Gaussian (with free cr and free centroid energy) 

narrow Gaussian and to 240 eV (the value in the second line in conclude that a double Gaussian is a better fit, but only at the 2cr 
Table 3) in the second broad Gaussian. level. 



We calculated the EW of the narrow and broad components, 
obtaining < 84 eV and 158*™ eV, respectively. This fit gives 
^ 2 /dof=19.15/26, to be compared with ^ 2 /dof =23. 70/27 of the 
fit with only narrow Gaussian (in Table 3). Hence, the signifi- 
cance of adding the broad component (Ay 2 =4.55 and Av=l) is 
> 95.4% but < 99%, calculated as mentioned in Sect. 4.2. We 



5. Conclusions 

We have studied the average spectrum of a sample of 123 
AGN with more than 200 net counts selected from the Chandra 
AEGIS, CDF-N and CDF-S surveys, covering < z < 3.5 and 
totalling >70000 cts. Compared to similar XMM-Newton and 
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Fig. 19. Fit of the average spectrum of the L — z subsamples (low L-low z subsamples in the left panels, low L- high z subsamples in 
the right bottom panel, high L-high z subsample in the right top panel) using the simulated continuum model + Gaussian (with free 
cr and free centroid energy) 
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Fig. 20. Dependence of unresolved 6.4 keV Fe line equivalent width with the average luminosity (left) and the average redshift (right) 
of the subsamples (first row of Table 3). Squares: L - z subsamples; circles: redshift subsamples; triangles: luminosity subsamples. 
We represent upper limits with a down-pointing arrow. Each set of points with the same symbol are statistically independent, but 
not so across symbols. 



Chandra studies of deep fields, we sample better the luminosity - 
redshift plane at z < 3 (including higher luminosities) and use 
deeper CDF-S data. With respect to shallower wider surveys us- 
ing XMM-Newton, we reach about an order of magnitude lower 
luminosities over the same redshift range. 

To improve the signal-to-noise ratio and to study the aver- 
age properties of the sample, we computed its average spectrum, 
adapting the method developed by |Corral et al. (2008)[ We con- 
structed the continuum model with simulations, in order to accu- 
rately analyze our resulting spectra. We also assessed the effect 



of our averaging method on unresolved features around 6.4 keV, 
to obtain the intrinsic widths of any putative lines. 

We repeated the averaging procedure for the subsamples de- 
fined in two intervals of column density, 2-10 keV X-ray lu- 
minosity, redshift, and in three 2D intervals in the luminosity- 
redshift (L - z) plane (see Table[T]l. 

We have estimated the significance of the presence of nar- 
row and broad features around 6.4 keV in a model-independent 
way (Table 2), finding that narrow features are significant at 
>98%, except for the high-L, high-z subsamples, where it is 
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lower. Broad features would be harder to detect and are most 
significant at low L and low z. 

Analyzing the spectra with XSPEC using the simulated con- 
tinuum and modelling the Fe line as a Gaussian, we have found 
strong evidence of a narrow Fe line at high significance in our 
full sample and in most of our subsamples, in particular (Table 
0: 

- We detected a strong (at > 99.73 %) narrow Fe line with 
EW= 85 + 35 eV and 73 + 32 eV in the absorbed and unab- 
sorbed subsamples, respectively (defined as having intrinsic 
column density above and below 10 21 5 cirT 2 , respectively), 
with a hint of higher central energy ~6.5 keV in the unab- 
solved case, but still compatible with 6.4 keV at 90%. 

- Segregating the sources purely on luminosity, we found that 
the average spectrum in the low-L (L < 8 x 10 43 erg s _1 ) 
subsample shows a significant and strong narrow Fe line 
(EW= 89+32 eV at >99.99%). The significance and strength 
of the high-L are clearly smaller: EW= 59 + 38 eV at 95.4%. 

- Separating instead at z = 1.005, the narrow Fe line is 
again stronger in the low-z subsample (EW= 80 + 32 eV 
at >99.99%) than in the high-z subsample (EW= 71 + 37 eV 
at >99.73%) but the difference is not significant 

- Defining distinct areas in the L - z plane (the division point 
was L - 14.3 x 10 43 erg s~',z = 0.76), the most significant 
line is found in the low L-low z subsample (EW= 108+42 eV 
at > 99.99%), with a hint of a broad component (<x = 240 + 
90 eV at > 2<x significance). The low L-high z subsample 
also shows a somewhat weaker narrow line (EW= 75 ±40 eV 
at ~ 99.7%). The Fe line detection is at < 90% significance 
in the high L-high z subsample (EW< 90 eV) 

We studied the trend of the EW for the narrow Fe line with 
redshift and luminosity of our subsamples. We did not find a sig- 
nificant dependence: in Fig.|20]that trend is not statistically more 
significant than a constant. In order to disentangle the effects of 
the two parameters on the line detection, surveys of AGN cover- 
ing a wider area in the L - z plane are needed. Since the Fe line 
is detected more significantly at low-L and low-z, samples of the 
local Universe are particularly suited to this end. 

We have also investigated whether more sophisticated mod- 
els for the line profile would provide a better fit, concentrating 
on the low L-low z subsample, since it is the only one where 
a broad line may have been detected. A diskline relativistic 
profile provides a worse fit than a broad Gaussian, with more 
parameters. Additionally, the inner radius of the disk would be 
too large to produce significant relativistic effects. Therefore, we 
do not find evidence for a relativistic profile in any of our average 
spectra. 

However, as has been pointed out in Guainazzi e t al. (2006)| 
in the average spectra the broad component of the observed line 
comes from a maximum ~ 50% of the sample: the line in an 
average spectrum cannot have a relativistic profile as definite 
as the one that can be in principle observed in a single, good 
quality spectrum. In particular, there can be a contribution from 
partially ionized Fe that produces narrow lines at energies > 6.4 
keV, up to ~1 keV if the Fe is completely ionized. Our observed 
~symmetric Gaussian profile may hence be a combination of 
many mildly ionized Fe lines with a red wing, perhaps from rel- 
ativistic effects. As mentioned in Sect. 1, a similar result has 
been recently found in the stacking of the X-ray spectra of the 
COSMOS sample dlwasawaet al. 201 11 1. 

Allowing for both a narrow (from neutral material far away 
from the central source, e.g. the putative torus) and a broad 



Gaussian (as just discussed) components provides a better fit 
than a single narrow or broad Gaussian, but with a modest sig- 
nificance < 99%. 

Further spectral averaging studies with higher statistics cov- 
ering a wider range of source properties would allow a more 
detailed characterization of the Fe feature in the X-ray spec- 
tra of AGN, and whether it depends on cosmic time, intrinsic 
brightness, amount of surrounding material, black hole mass, 
accretion rate, and other physical parameters, helping to con- 
strain the physical properties of AGN throughout the history of 
the Universe. 
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Appendix A: Results with different definitions of 
the sample 

We assessed the robustness of our results by changing the signal- 
to-noise ratio that characterizes our sample. To do that, we re- 
peated the analysis for all the sources with more than 50 and 
100 counts in 2 - 12 keV. We excluded from these samples the 
spectra with the lowest continuum flux between 2-5 keV (a bias 
during the normalization process can be introduced, see Sect.3.1 
for the details for method): CDFS-227 (CDF-S source with RA: 
53.082, DEC: -27.690), CDFN_405 (CDF-N source with RA: 
189.431, DEC: 62.177), EGS1-003 (AEGIS source with RA: 
215.76, DEC: 53.45). After having defined the samples, we re- 
defined the subsamples of intrinsic column density, luminosity 
and redshift, as described in Sect. 2.3. The results are in Table 
I A. II and [A. 21 We can see from the last column in the tables that 
the Fe line significance grows with the signal-to-noise ratio of 
the sample. After having applied the same analysis method used 
for the default > 200 cts sample, we repeated the analysis in 
XSPEC for the > 50 cts and > 100 cts sample, finding consistent 
results. 

We have made the simulations of an unresolved 6.4 keV line 
for the > 50 cts sample, obtaining results compatible with that 
of the default > 200 cts sample. 

Appendix B: Tests of the method: approximation 
made in the correction of the spectra for 
detector response 

We checked whether the correction for the response matrices de- 
pends strongly on the model used to do it. To make this check, 
we ran XSPEC simulations. 

We start with an "average" but high signal-to-noise ratio 
source modelled with an intrinsically absorbed powerlaw and 
an unresolved Gaussian line at 6.4 keV, with the average val- 
ues of the full sample (see Table Q]): z = 1.15, T = 1.24, 
N H = 1.75 x 10 22 ctrr 2 , EW=100 eV. We compared the input 
model with the response-matrix-corrected spectra ("unfolded" 
spectra in XSPEC) assuming three different models: 

1 . the input model without the Gaussian 

2. a powerlaw with F = 

3. a powerlaw with F = 2. 

We can see the results in Figs. lB. ll in the observed frame. The 
best approximation is obtained using the input parameters, as ex- 
pected. However, in a wide range around the line, the result does 
not depend noticeably on the model used for the response ma- 
trix correction, even for models clearly different from the input 
one. The significant discrepancies start only below about 2 keV 
rest-frame. 

We ran another simulation for an absorbed source, using the 
continuum parameters of one source in the absorbed subsample, 
(z = 1.15, F = 1.64; N H = 43.22 x 10 22 cirT 2 , ) and a Gaussian 
at 6.4 keV with cr — 0, EW = 100 eV, again with high signal- 
to-noise ratio. We show in Fig. lB.2l the input model and the "un- 
folded" spectra using the three models above. In this case we 
also see that there is no important difference in the continuum 
close to the Fe line between the various models. The differences 
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in ^ ^ ^ ^ 

d0.5 1 2 5 

Energy (keV) 



Fig. B.l. Comparison, for a simulated "average" source, be- 
tween unfolded spectrum with best fit (green circles), with 
gamma =2 (red plus), with gamma =0 (continuous black line), 
with model (dashed blue line). 




^0.5 1 2 5 



Energy (keV) 

Fig. B.2. Comparison, for a simulated absorbed source, between: 
unfolded spectrum with best fit (green circles), with gamma =2 
(red plus), with gamma =0 (continuous black line), with model 
(dashed blue line). The unfolded spectrum obtained with the best 
fit is along the model, with the exception of the Fe line region, 
where it deviates. 





* 





























2 5 10 



Energy (keV) 

Fig. B.3. Comparison between the average spectra after having 
corrected for the detector response using: the best fit (average 
observed spectrum represented by full circles, simulated spec- 
trum by a dashed line) and with the powerlaw with Y - 2 (aver- 
age observed spectrum represented by stars, simulated spectrum 
by a contiuous line). Sample with the sources with more than 
200 net counts. 

for the "unfolding" appears to produce a more conservative esti- 
mate of the Fe line flux. 



below ~3 keV rest-frame are much stronger in this case com- 
pared to the "average" source, but again the result of the input 
model follows very closely the input spectrum. 

Therefore, we chose to correct the spectra for the response 
matrices using their own best-fit continuum models and we did 
not simply divide the spectra for the effective area of the detector. 
This simple method, actually assumes a flat continuum and does 
not account for the limited spectral resolution of X-ray detectors. 

It is important to notice that in any case there is an energy 
blending that widens the Fe line. We have quantified this effect 
using simulations, as explained in Sect. 3.3. 

In order to check the effect of the correction for detector re- 
sponse on the average spectrum, we corrected the spectra for 
the detector responses using a powerlaw with F = 2. We then 
made the average of the observed and simulated spectra, using 
the same procedure explained in the Sect. [3] We can see in Fig. 
IB.3l the results: the differences are again small and, in any case, 
within the mutual error bars. If anything, using the best-fit model 
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Table 1. Properties of the full sample and the subsamples (see text). 



Sample 


N 


#2-12 


#5-8 


<z> 


<^43> 


W,22> 


r 


* 2 /dof 


r 

A sim 


@~ sim 


EW S ^, 












10 4j erg s 1 


10 22 cm 2 








eV 


eV 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


CDF-S 


33 


21561 


5223 


1.26 


5.20 


3.80 


- 


- 


- 


- 


- 


CDF-N 


25 


16134 


3429 


1.11 


8.00 


0.88 


- 


- 


- 


- 


- 


AEGIS 


65 


32971 


6873 


1.11 


20.31 


1.04 


- 


- 


- 


- 


- 


Total 


123 


70667 


15526 


1.15 


13.82 


1.75 


1.24 + 0.05 


9.77/13 


1.897 ±0.001 


117 ± 1 


182 


log(N H ) > 21.5 


54 


27888 


7304 


1.15 


14.11 


3.93 


0.48 ± 0.07 


8.94/13 


1.898 ±0.001 


129 ± 1 


190 


log(N H ) < 21.5 


69 


42778 


8222 


1.14 


13.55 


0.05 


1.78 + 0.06 


14.96/13 


1.896 ±0.001 


116 ± 1 


181 


L 43 < 8 


74 


34620 


7630 


0.85 


2.89 


1.64 


1.21 ± 0.06 


10.42/13 


1.893 ± 0.001 


111 ± 1 


183 


L 43 > 8 


49 


36047 


7896 


1.60 


30.20 


l.9l 


1.31 ±0.07 


14.70/13 


1.878 ±0.001 


128 ± 1 


174 


z < 1.005 


63 


34843 


7507 


0.70 


5.43 


1.28 


1.23 ±0.06 


10.57/13 


1.896 ±0.001 


105 ± 1 


186 


z > 1.005 


60 


35824 


8019 


1.62 


22.50 


2.24 


1.24 ±0.07 


5.66/13 


1.886 ±0.001 


138 ± 1 


180 


z < 0.76 


39 


23713 


5285 


0.57 


3.23 


1.17 


1.31 ±0.07 


11.68/13 








L 43 < 14.3, z < 0.76 


36 


19007 


4446 


0.57 


2.00 


1.25 


1.29 ±0.08 


10.99/13 








L 43 < 14.3, z > 0.76 


49 


24286 


4902 


1.14 


5.58 


1.90 


1.11 ±0.08 


7.27/13 








L 43 > 14.3, z > 0.76 


35 


22669 


5339 


1.80 


37.05 


2.19 


1.34 ±0.09 


11.76/13 









Notes. Columns: (1): (Sub)Sample; (2): Number of sources; (3): number of net counts in 2-12 keV rest-frame; (4): number of net counts in 5-8 
keV rest frame; (5): average redshift; (6): average rest-frame 2-10 keV luminosity in units of 10 43 erg s~', corrected for Galactic and intrinsic 
absorption*; (7): average intrinsic column density in 10 22 crrT 2 *; (8): gamma from the fit of the average observed spectrum with a powerlaw 
between 2 and 5 keV (rest-frame); (9):^ 2 /dof of the fit of the average observed spectrum with a powerlaw between 2 and 5 keV (rest-frame); (10): 
slope of the powerlaw ** of the average spectrum obtained with the simulations of the Fe line (see text); (11): width obtained in those simulations; 
(12): EW of the line obtained from those simulations. * Obtained averaging the values obtained from the fit to the individual spectra (see text); ** 
obtained fitting it with a powerlaw plus a Gaussian Fe line 
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Table 3. Results of fits of the average spectrum of the full sample and its subsamples. 



Sample 


* 2 Mof(s) 


* 2 /dof(g) 


Pi 


Pi 


E 


cr 


EW 


(z) 




N h ,22 








% 


% 


keV 


eV 


eV 




l0 4J erg s 1 


\0 12 cm 1 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


Total 


47.50/28 


23.74/27 


> 99.99 




6.4 





74 ± 24 


1.15 


13.80 


1.75 






21.98/26 


> 99.99 


< 90 


6.4 


< 138 


94 ± 30 












23.74/26 


> 99.99 


< 90 


6.40 ±0.10 





73 ± 25 












21.98/25 


> 99.99 


< 90 


6.40 ± 0.10 


< 123 


92 ± 28 








1 /XT \ . "11 — 

log(N H ) > 21.5 


35.61/28 


21.66/27 


> 99.73 




6.4 





85 ± 35 


1.15 


14.11 


3.93 






21.15/26 


> 99.73 


< 90 


6.4 


< 123 


100 ± 45 












21.42/26 


> 99.73 


< 90 


6.28 ± 0.02 





92 ± 41 












21.22/25 


> 99.73 


< 90 


6.40 + 0.30 


< 151 


91 ± 39 








log(N H ) < 21. 5 


53.03/28 


39.38/27 


> 99.73 




6.4 





73 ± 32 


1.14 


13.55 


0.045 






38.18/26 


> 99.73 


< 90 


6.4 


< 189 


100 ± 43 












37.41/26 


> 99.73 


< 90 


6.50 ±0.10 





85 ± 36 












37.15/25 


> 99.73 


< 90 


6.50 ±0.10 


< 59.6 


87 ± 36 








L 43 < 8 


44.26/28 


23.96/27 


> 99.99 




6.4 





89 ± 32 


0.85 


2.89 


1.64 






23.70/26 


> 99.99 


< 90 


6.4 


< 72 


97 ± 35 












23.96/26 


> 99.99 


< 90 


6.4 ±0.1 





89 ± 32 












23.56/25 


> 99.73 


< 90 


6.4 ± 0.9 


< 68 


97 ± 34 








L 43 > 8 


42.59/28 


36.42/27 


> 95.4 




6.4 





59 ± 38 


1.597 


30.20 


1.91 






36.41/26 


~ 95.4 


< 90 


6.4 


< 7 


58 ± 39 












36.42/26 


> 90.0 


< 90 


6.4 ± 0.2 





58 ± 39 












34.86/25 


> 90.0 


< 90 


6.30 ± 0.20 


< 7 


75 ±45 








z < 1.005 


52.34/28 


35.26/27 


> 99.99 




6.4 





80 ± 32 


0.704 


5.43 


1.28 






34.21/26 


> 99.73 


< 90 


6.4 


< 100 


94 ± 37 












33.69/26 


> 99.99 


< 90 


6.50 ± 0.10 





92 ± 35 












33.28/25 


> 99.73 


< 90 


6.45 ± 0.10 


< 101 


99 ± 37 








z > 1.005 


29.56/28 


19.77/27 


> 99.73 




6.4 





71 ± 37 


1.617 


22.5 


2.24 






19.45/26 


> 99.73 


< 90 


6.4 


< 120 


86 ± 44 












19.77/26 


> 99.73 


< 90 


6.40 ±0.15 





71 ± 37 












18.66/25 


> 99.73 


< 90 


6.30 ± 0.10 


< 76 


85 ± 43 








L 43 < 14.3, z < 0.76 


42.01/28 


23.70/27 


> 99.99 




6.4 





108 ± 42 


0.566 


2. 


1.254 






19.23/26 


> 99.99 


> 95.4 


6.4 


240 ± 90 


178 ± 63 












23.70/26 


> 99.73 


< 90 


6.40 ± 0.15 





108 ± 42 












19.07/25 


> 99.99 


> 90 


6.42 ± 0.12 


216 ± 50 


172 ± 60 








z < 0.76 


43.28/28 


23.04/27 


> 99.99 




6.4 





108 ± 40 


0.573 


3.23 


1.173 






19.04/26 


> 99.99 


~ 95.4 


6.4 


157 ± 30 


154 ± 50 












23.04/26 


> 99.99 


< 90 


6.40 ± 0.15 





109 ± 39 












18. /2/z5 


> 99.99 


< 90 


6.43 ±0.11 


215 ± 50 


170 ± 56 








T 1/11 -T ■ — C\ H£* 

L43 < 14.3, z > U. 10 


nnno 

zj.y //zo 


lo.yy/z/ 


no n 
~ yy. 1 




0.4 


U 


lj± 4U 


1 i/ii 
1. 141 


3. 38 


1 onn 






16.99/26 


> 95.4 


< 90 


6.4 


< 4 


75 ±40 












16.99/26 


> 95.4 


< 90 


6.40 ±0.10 





75 ±40 












16.83/25 


>95.4 


< 90 


6.43 ±0.10 


< 5 


76 ±42 








L43 > 14.3, z > 0.76 


31.85/28 


29.56/27 


<90.0 




6.4 





< 90 


1.803 


37.05 


2.187 






29.56/26 


< 90.0 


< 90 


6.4 


< 5 


< 91 












29.09/26 


< 90.0 


< 90 


6 05 +0 ' 60 

u - u -0.05 





50 ±47 












26.72/25 


< 90.0 


< 90 


6.22 ±0.15 


<5 


73 ±53 









Notes. L43: luminosity in units of 10 43 erg s -1 . Columns: (1): sample; (2): ^ 2 /dof of the fit with simulated continuum model; (3): ^ 2 /dof of the 
fit with the same continuum and the Gaussian; (4): probability P(Ax 2 , Av)s, g (of adding the Gaussian to the model, see text); (5): probability 
P(Ax 2 , Av)g, go (of leaving the parameter of the Gaussian free, see text); (6): central energy of the Gaussian; (7): a after subtracting the intrinsic 
sigma obtained in the simulations of the Fe line; (8): EW of the Gaussian; (9): average redshift of the sample; (10): average luminosity of the 
sample in (10 43 ) erg s" 1 ; (11): average column density of the sample in 10 22 cm' 2 . 



16 



S. Falocco et al.: Averaging the AGN X-ray spectra from deep Chandra fields 
Table A.l. Properties of the sample composed by all the sources with more than 50 counts. 



Sample 


N 


N2-12 


Ns-i 


(z) 


<^43> 

10 43 erg s _1 


10 22 x cnr 2 


Significance 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Full 


347 


93385 


21730 


1.20 


8.60 


3.08 


85 


log(N H ) > 21.5 


177 


40174 


10944 


1.14 


8.00 


6.04 


85 


log(N H ) < 21.5 


170 


53211 


10786 


1.28 


9.50 


0.04 


72 


L 43 < 6 


234 


46546 


11021 


0.90 


1.60 


3.09 


89 


L 43 > 6 


113 


46840 


10709 


1.80 


23.40 


3.07 


74 


z < 1 


184 


46672 


10530 


0.65 


2.60 


3.08 


83 


z > 1 


163 


46713 


11199 


1.83 


15.60 


3.08 


85 


z < 0.76 


120 


31537 


7511 


0.53 


1.80 


3.75 


81 


L43 < 8.4, z < 0.76 


116 


26743 


6614 


0.52 


1.20 


3.35 


82 


L43 < 8.4, z > 0.76 


157 


31862 


6866 


1.34 


1.70 


2.47 


94 


L43 > 8.4, z > 0.76 


70 


29985 


7353 


2.08 


32.10 


3.32 


77 



Notes. L 43 : luminosity in units of 10 43 erg s -1 . Columns: (1): (Sub)Sample; (2): Number of sources; (3): Number of counts in 2 -12 keV; (4): 
Number of counts in 5-8 keV; (5): average redshift; (6): average luminosity in 10 43 erg s -1 ; (7): average column density of the local absorber in 
10 22 - cm" 2 ; (8): significance of the Fe line calculated as the number of average simulated spectra with a lower flux than the flux of the average 
observed spectrum (calculation made between between 6.2 and 6.6 keV). 

Table A.2. Properties of the sample composed by all the sources with more than 100 counts. 



Sample 


N 


#2-12 




(z) 


(U 3 ) 




Significance 












10 43 erg s" 1 


10 2 W 2 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Full 


219 


84275 


19312 


1.24 


11.83 


2.63 


92 


log(N H ) > 21.5 


104 


34893 


9422 


1.15 


10.80 


5.50 


91 


log(N H ) < 21.5 


115 


49382 


9890 


1.32 


12.69 


0.04 


100 


L43 < 6.65 


134 


42061 


9718 


0.91 


2.25 


3.03 


91 


L43 > 6.65 


85 


42214 


9594 


1.76 


26.82 


1.99 


94 


z < 1.01 


109 


41477 


9216 


0.68 


3.80 


2.74 


90 


z > 1.01 


110 


42798 


10096 


1.80 


19.71 


2.52 


99 


L 43 < 21.0, z< 0.78 


71 


28349 


6587 


0.55 


2.40 


3.58 


90 


L43 <21.0, z>0.78 


93 


29478 


6273 


1.30 


4.80 


2.16 


100 


L 43 > 21.0, z > 0.78 


55 


26448 


6452 


2.02 


35.80 


2.20 


90 



Notes. L 43 : luminosity in units of 10 43 erg s -1 . Columns: (1): (Sub)Sample; (2): Number of sources; (3): Number of counts in 2 -12 keV; (4): 
Number of counts in 5-8 keV; (5): average redshift; (6): average luminosity in 10 43 erg s -1 ; (7): average column density of the local absorber in 
10 22 - cirr 2 ; (8): significance of the Fe line calculated as the number of average simulated spectra with a lower flux than the flux of the average 
observed spectrum (calculation made between between 6.2 and 6.6 keV). 
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